BCS-BEC quantum phase transition and collective excitations in two-dimensional Fermi gases 

with p- and J-wave pairings 
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It is generally believed that the BCS-BEC evolution in fermionic systems with s-wave pairing is a smooth 
crossover. However, for nonzero orbital-angular-momentum pairing such as p- or d-wave pairing, the system 
undergoes a quantum phase transition at the point where the chemical potential p vanishes. In this paper, we 
study the BCS-BEC quantum phase transition and the collective excitations associated with the order-parameter 
fluctuations in two-dimensional fermionic systems with p- and rf-wave pairings. We show that the quantum 
phase transition in such systems can be genetically traced back to the infrared behavior of the fermionic ex- 
citation at p = 0: £k ~ k 1 , where / = 1,2 is the quantum number of the orbital angular momentum. The 
nonanalyticity of the thermodynamic quantities is due to the infrared divergence caused by the fermionic ex- 
citation at p = 0. As a result, the evolution of the Anderson-Bogoliubov mode is not smooth: Its velocity is 
nonanalytical across the quantum phase transition. 

PACS numbers: 03.75.Ss, 05.30.Fk, 67.85.Lm, 74.20.Fg 



I. INTRODUCTION 

It was proposed by Eagles Utl and Leggett several 
decades ago that, by tuning the strength of the attractive in- 
teraction in a many-fermion system, one can realize an evolu- 
tion from the Bardeen-Cooper-Schrieffer (BCS) superfluidity 
to Bose-Einstein condensation (BEC) of difermion molecules. 
For s-wave interaction, it is generally believed that the BCS- 
BEC evolution is a smooth crossover ll3l- fl0ll . One of the most 
interesting systems is the dilute spin-| Fermi gas in three di- 
mensions with short-range s-wave attraction. The system is 
characterized by a dimensionless parameter 1 / (kpa s ), where 
a s is the s-wave scattering length of the short-range interaction 
and kp is the Fermi momentum in the absence of interaction. 
For such a system, one finds a smooth BCS-BEC crossover 
when the parameter l/(k F a s ) goes from -co to oo. In addi- 
tion, the Anderson-Bogoliubov collective mode of fermionic 
superfluidity with weak attraction evolves smoothly to the Bo- 
goliubov excitation of weakly repulsive Bose condensate with 
strong attraction JH [HI O . The smooth BCS-BEC crossover 
with s-wave interaction has been experimentally studied by 
using ultracold fermionic atoms lfl3ll . where the s-wave scat- 
tering length was tuned by means of the Feshbach resonance. 

On the other hand, for nonzero orbital-angular-momentum 
pairing, such as p- or of- wave pairing Ill4l - tl9ll . it is generally 
accepted that the BCS-BEC evolution is not smooth but as- 
sociated with some quantum phase transition. Such quantum 
phase transition cannot be characterized by a change of sym- 
metry or the associated order parameter. Instead, different 
quantum phases can be distinguished topologically fTill . The 
quantum phase transitions and the finite temperature phase di- 
agrams of single-species polarized Fermi gases tuned across 
a p-wave Feshbach resonance were studied by Gurarie, Radz- 



ihovsky, and Andreev IU7H (for an extensive study, see Ref. 
ifioll ). In three dimensions, they showed that across the p- 
wave Feshbach resonance, the system would undergo a quan- 
tum phase transition from a p^-wave to a p x + ;/?,,-wave su- 
perfluid. Furthermore, the latter state undergoes a topological 
transition at zero chemical potential p. 

In this paper, we explore some generic aspects of the BCS- 
BEC quantum phase transition in two-dimensional (2D) p- 
and of-wave paired Fermi superfluids, which is expected to oc- 
cur at the point where the chemical potential p vanishes ifTil - 
[l6ll . While for nonzero orbital-angular-momentum pairing it 
is important to explore the order parameter symmetry, in thi s 
paper we show that the generic nature of the BCS-BEC quan- 
tum phase transition at p — in two-dimensional fermionic 
systems is independent of the order parameter symmetry but 
depends on the infrared behavior of the order parameter. 

In general, the superfluid order parameter or the gap 
function A(k) is momentum dependent for nonzero orbital- 
angular-momentum pairing. The crucial observation in this 
paper is that the infrared behavior of the order parameter A(k) 
depends only on the quantum number / of the orbital angular 
momentum. For /th-wave pairing, we have 



A(k) ~ ie 



0. 



(1) 
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For instance, the order parameter of p-wave pairing can be 
expressed as A(k) = Aokg(tp) where <p is the polar angle in 
two-dimensions. The anisotropic p x (p v ) pairing corresponds 
to g(ip) = cos ip (sin^) and the complex p x + ip y pairing cor- 
responds to g(tp) = For any case, we find that the infrared 
behavior of the order parameter is A(k) ~ k, independent of 
the order parameter symmetry. 

The single-particle excitation spectrum in the superfluid 
state reads E k = [£ 2 +|A(k)| 2 ] 1/2 , where £ k = k 2 /(2M)-p, with 
M being the fermion mass. Therefore, for p-wave (/ = 1) and 
cZ-wave (/ = 2) pairings, the infrared behavior of the single- 
particle spectrum at the quantum phase transition point p = 



2 



is solely determined by the order parameter, that is 
EkO" = 0) ~ kf, k -> 0. 



(2) 



Note that the above behavior applies only to I — 1 and / = 2. 
For / > 2, the infrared behavior becomes dominated by the 
kinetic term In this paper, we show that the type of 2D 
momentum integral (m, n integer) 



Jr-oo in 
kdk — 



(3) 



appears in the expressions of various physical quantities. This 
integral is divergent in infrared at fi = if Im > n + 2. 
Therefore, for a 2D system with nonzero orbital-angular- 
momentum pairing, nonanalyticities generally appear due to 
the infrared divergence at the quantum critical point jj — 0. 

The nonanalyticities of the thermodynamic quantities at 
vanishing chemical potential which signal the quantum phase 
transitions in 2D p- and d-wave paired superfluids were first 
studied by Botelho and Sa de Melo using Nozieres-Schmitt- 
Rink (NSR) type potentials lfl5l, [l6ll . In this paper, we extend 
their conclusions significantly. We show that the nonanalyt- 
icities are essentially related to the infrared behavior of the 
interaction potential. Therefore, the quantum phase transition 
at vanishing chemical potential is independent of the details of 
the interaction potential as well as the order parameter sym- 
metry, which can be intuitively understood by the integral <j3j . 
We also show that the evolution of the Anderson-Bogoliubov 
mode is also not smooth; that is, the sound velocity goes non- 
analytically across the quantum phase transition. 

In the rest of this paper, we study the BCS-BEC quantum 
phase transition and the behavior of Anderson-Bogoliubov 
mode across the phase transition for two typical cases: p-wave 
pairing in spinless Fermi gases in Sec. II and cZ-wave pairing 
in spin-i Fermi gases in Sec. III. 



II. p-WAVE PAIRING IN SPINLESS FERMI GASES 

For a 2D interaction potential V(r), the momentum-space 
matrix element V(k, k') can be expressed as lfl5ll 



V(k, k') = ^ e'" flkk ' V (n, (k, k'), 



(4) 



where #kk' is the angle between k and k' and the coefficients 
V^ik, k') are given by 



V (n \k,k') =2tt drrJ n {kf)J n (k'r)V(r). 



(5) 



Here J„(x) is the Bessel function of order n. The « = +1 
components correspond to the Ith angular momentum channel. 

It is possible to retain only the n = +1 terms, by assum- 
ing that the dominant contribution to the scattering processes 
between fermions occurs in the Ith angular momentum chan- 
nel. This assumption may be experimentally relevant since 
atom-atom dipole interactions split different angular momen- 
tum channels such that they may be tuned independently. 



Using the properties of the Bessel function, in the low en- 
ergy limit k — > and k' — > we have 



V (±t \k,k')~k'(k'y. 



(6) 



Therefore, the potential becomes separable in the low energy 
limit. While it is certainly not separable for general values of 
k and k', in the following we use separable potential for the 
Zth wave interaction to simplify our formulation. Since the 
generic features of the BCS-BEC quantum phase transition 
with nonzero orbital-angular-momentum pairing are related 
solely to the infrared behavior of the interaction potential, the 
use of a separable potential is without loss of generality. 

The many-body Hamiltonian of 2D spinless fermions can 
be written as H = Ho + H mt , where the single-particle part 
reads Hq - £k £k<A k <A k and the p-wave interaction part can be 
written as 



flint = J] ^( k ' k ')'AL q /2 i / r -k +q /2 ( /'-k' +q /2 i / r k' + q 



(7) 



k.k'q 



The generic infrared behavior of the p-wave interaction poten- 
tial is V p (k, k') ~ kk' for k, k' — » 0. Without loss of generality, 
we consider a separable potential for the p-wave interaction, 
V p (k, k') = -/tpr p (k)r*(k'), where A p is the coupling con- 
stant and the function T p (k) characterizes the p-wave pairing 
symmetry. The infrared behavior of the function F ; ,(k) is 



T p (k) ~ k, k -> 0. 



(8) 



In the functional path integral formalism, the partition func- 
tion of the system at finite temperature T is given by Z, = 
j D\pD\ji^ exp {-J>[^, (^]J, where the action reads 



*(TMT) + fl(^^) 



(9) 



To decouple the interaction term we introduce an auxil- 
iary complex pairing field <£> q (T), which couples to (// i/z^ 
and apply the Hubbard-Stratonovich transformation. Us- 
ing the Nambu-Gor'kov basis ^ = (>^ k ,i^l) T and integrat- 
ing out the fermionic degrees of freedom, we obtain X = 
exp { - ^ e ff[ < l ) , < I ) *]), with the effective action given 

by 



5eff = r dr^^ft^ " TrlnG- ) 



(10) 



where the inverse single-particle Green's function G kk , 
given by 



is 



G k,k' 



-Or + &)<5k,k< %-k'(T)r p (^) 



-(dr-£k)5k,k< 



(ID 



The effective action S e s [Q>, O*] cannot be evaluated pre- 
cisely. In this work, we consider mainly the zero-temperature 
case. Therefore, we follow the conventional approach to the 
BCS-BEC crossover problem; that is, we first consider the su- 
perfluid ground state which corresponds to the saddle point 



3 



of the effective action and then study the Gaussian fluctuation 
around the saddle point. The Gaussian-fluctuation part corre- 
sponds to the collective modes, including the gapless Gold- 
stone mode or the Anderson-Bogoliubov mode. This Gold- 
stone mode appears as a resonance in the spectrum of the 
density-density correlation function lEoll . Therefore, it can 
be experimentally probed by many techniques, including two- 
photon Bragg scattering 1I21I l22ll . 

In the superfluid ground state, the pairing field <1> acquires 
a nonzero expectation value Ao, which can be set to be real 
without loss of generality. Then we decompose the pairing 
field as O q (r) = Ao5 q ,o + (pqir), where (f> is the fluctuation 
around the mean field. The effective action S e ff [<5, <5*] can be 
expanded in powers of the fluctuation 0, that is, 



>S eff [0,0*] 



*] + 



(12) 



where 5S(Ao) is the saddle point or mean-field effective 
action with Ao determined by the saddle point condition 

as^/a&o = 0. 



A. Quantum phase transition 

Neglecting the beyond-mean-field contribution, which is 
generally thought to be small, we obtain the effective potential 



^(A ,ju) 



E k ) 



(13) 



at zero temperature, where the single-particle excitation spec- 
trum Ek is given by 



£ k = yjil + |A(k)| 2 . 



(14) 



Here A(k) = AoTy^k) is the gap function for /j-wave pairing. 
The order parameter Ao is determined by minimizing the ef- 
fective potential, which gives the gap equation 



ir,(k)p 

4£ k 



(15) 



To study the nature of the quantum phase transition, we 
focus on the thermodynamic potential Q(,u) = r V(Ao(fi),/i). 
Note that the thermodynamic variable here is the chemical po- 
tential fi. The order parameter Ao should be solved as an im- 
plicit function of ft through the gap equation. The first deriva- 
tive of the thermodynamic potential with respect to fx gives the 
number equation, 



n(ji, A ) 



<9Q(yU) 

dfi 



&('-!)■ 



Then we consider the second derivative 

d 2 Q.(fi) 



dfi 2 



(16) 



(17) 



Considering the fact that Ao is an implicit function of ft, we 
obtain 



dn(n,A ) dn(n,Ao)dA 
a(ji) = + 



dji 



8A Q dfi 



(18) 



The derivative of Ao with respect to fi can be obtained from 
the gap equation. We have 



3Ao _ dn( M ,A )(d 2 'V(Ao,n) 
dfi 



dAo 



dAl 



Therefore, a(/j) can be expressed as 



atju) 



dn(jx, A ) / dn(pi, A ) 



dfi 



dAo 



dAl 



(19) 



(20) 



We note that the second term was missing in the previous 
study lfl5h . However, the nonanalytical behavior is dominated 
by the first term. Therefore, the conclusions in Ref. lfi"5ll are 
still reliable. The derivatives in Eq. (20) can be explicitly 
evaluated as 



dnjfi, A ) _Ajy 
dfi = 2 2u 

dn(jx, A ) A 



dAo 



- z2 V 

" 2 V 



ir P (k)| 2 

E 3 ' 

^kir„(k)| 2 



d 2 V{A ,fi) A 2 



OA 2 



^k 

T P (k)| 4 

E 3 ' 
c k 



(21) 



By using an NSR-type potential, Botelho and Sa de Melo 
found that the quantity a (proportional to the isothermal com- 
pressibility) is nonanalytical at fi — Here we show 
that this nonanalyticity generally appears due to infrared di- 
vergence at fi — 0. Actually, for k — > 0, we have F ; ,(k) ~ k 
and Ek ~ k at fi — 0. Therefore the momentum integrals in 
Eq. (18) are infrared safe at fi = 0. Then we further consider 
the derivatives of a(ji) with respect to fi. Actually, without 
explicit calculations, we find that the following momentum 
integral 



T P (k)p 

E 5 



(22) 



appears in the expression of d 2 a{fi)/dfi 2 . At vanishing chem- 
ical potential fi = 0, the infrared behavior of the integral is 



Ji(0) 



r kd k-~ r 

Jo k 5 Jo 



dk 



(23) 



Therefore, this integral is infrared divergent. As a result, the 
derivative d 2 a(fi)/dfi 2 is divergent at the point fi = 0, which 
indicates that da(fi)/dfi is discontinuous at fi = 0. Thus the 
quantity a(fi) itself is nonanalytical across the quantum criti- 
cal point fi — 0. 

The above discussion shows that a quantum phase transi- 
tion occurs at fi — 0, indicated by the nonanalyticity of the 
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thermodynamic potential. The nonanalyticity is caused solely 
by the infrared behavior of the p-wave pairing potential and 
is independent of the details of the interaction as well as the 
order parameter symmetry (i.e., the angle dependence of the 
gap function). Therefore, the BCS-BEC quantum phase tran- 
sition is quite generic in two-dimensional systems, driven by 
the infrared divergence at p — 0. 

In Fig. Q] we show the numerical results of a(p) for two 
typical p-wave pairings: the complex and isotropic p x + ip v 
pairing and the anisotropic p x pairing. For the p x +ip y pairing, 
the superfluid state is fully gapped for both p > and p < 0, 
and it is only gapless at p — 0. For p x pairing, the superfluid 
state is gapped for p < and gapless for p > 0. In the cal- 
culation we have adopted a NSR-type potential to regularize 
the ultraviolet divergence in the gap equation. However, as we 
have argued above, the nonanalytical behavior at p = does 
not depend on this specific choice. The results in Fig. [T]show 
that the nonanalytical behavior appears for both isotropic and 
anisotropic pairings, consistent with our conclusion that the 
nonanalytical behavior is solely due to the infrared behavior 
of the pairing interaction. 



B. Collective modes 

To study the behavior of the collective modes across 
the BCS-BEC quantum phase transition, we consider the 

(21 

Gaussian-fluctuation part 5 <p*]. It can be written in the 
bilinear form 

^ff = \ X At (e)M(0A(0, (24) 
Q 

where Q = (iv„,q) with v„ being the boson Matsubara fre- 
quency, A(0 = [0(0, 0*(-0] T , and the 2 x 2 matrix M(0 
is the inverse of the collective-mode propagator. The matrix 
elements of M(0 are constructed by using the mean-field 
fermion propagator @{K) which is obtained from the expres- 
sion of Gk.k'M- Here K = {ico„, k), with a>„ being the fermion 
Matsubara frequency. We have 

Then the matrix elements of M(0 can be expressed as 

Mn(0 = ■i + iy'|r„(k)| 2 Tr[^ 11 (^ + W 2 2(^-)], 
M 12 (0 = ^|F^k)| 2 Tr[£ 12 (tf + )£ 12 (^)], (26) 

K 

where K ± - K + Q/2. In addition, we have M 2 2(0 = 
Mn(-0 and M 2 i(0 = M i2 (0. Completing the summa- 
tion over the boson Matsubara frequency, we obtain at zero 




9 ^ 
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0.5 1 1.5 

n 

FIG. 1 : The second derivative of the thermodynamic potential Q(p) 
with respect to p as a function of the density n for the isotropic p x +ip y 
pairing (upper panel) and the anisotropic p x pairing (lower panel). 
The inserts show the chemical potential p as a function of n. In this 
calculation, we employ the NSR-type potential 0] for the 2D case, 
T(k) = (klh)g( V )l(\ + k/k ) 3 ' 2 CI. Here g{tp) = e"" for p, + ip y 
pairing and g(tp) = cos<p for p x pairing. In this potential model, 
£o ~ r^ 1 , where ro plays the role of the interaction range, ko and 
ki set the momentum scales in the short and long wavelength limits, 
respectively. The density n has been normalized such that n = 1 at 
p = 0. All quantities are shown in proper units. 



temperature 
M„(0 

= 1 , y Ilp(k)| 2 / u\ul v 2 y_ \ 
M 12 (0 

Z |r p (k)| 2 / u+v+u-v- u + v + u-v- \ 

2 \iv„ -£+-£■_ iv„ +£■++£'_/' 

Here + and - denote the momenta k+q/2 and k-q/2, respec- 
tively, and are the BCS distribution functions defined as 

= = (1/2)(1 + &/£*)■ 

Taking the analytical continuation iv n — > a> + i0 + , the 
dispersions u>(q) of the collective modes are determined by 
the equation detM[w(q), q] = 0. It is usual to decompose 
M n (u),q) as Mn(w,q) = M+(w,q) + M~(w,q), where 
MJ"j (a), q) and M n (w, q) are even and odd functions of <d, 
respectively. Meanwhile Mi 2 (w, q) and M 2 i(tL>, q) are even 
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functions of a> automatically. Then we separate the complex 
field 4>{Q) into its amplitude part A(Q) and phase part 9(Q), 
<t>(Q) = MQ) + i'Ao0(0. The effective action 5® becomes 

^4^(^(0 <T(fi))N(fi)(jg), (28) 

where the matrix elements of N(0 read Nn(0 = 2(M|j + 
Mia), N 22 (0 = 2Ag(M+ - M 12 ), N 12 (0 = 2iA M" , and 
N 2 i(0 = -2;A M n . Since M n (0,q) = 0, the amplitude and 
phase modes decouple completely at a> = 0. At the saddle 
point we have precisely M|j(0, 0) = Mi 2 (0, 0). Therefore, the 
phase mode at q = is gapless, that is, the Goldstone mode 
or the Anderson-Bogoliubov mode for fermionic superfluids. 

To study the low-energy behavior of the collective modes, 
we make a small q and lo expansion of N(0 at zero tem- 
perature. In general, the expansion takes the form Nn = 
A + Cq 2 - Dto 2 + ■ • • , N 22 = Jq 2 - Rcu 2 + ■ • • , and N i2 = 
-N 2 i = -iBto + ■ ■ ■ . The explicit forms of the expansion pa- 
rameters A, B, D, and R can be evaluated as 



Therefore, we have 



4) y 
2 V 



ir P (k)| 4 



B = 



D = 



Ap y &ir p (k)p 

4 V 



£ 3 



i y ffir p (k)i 2 

8 ^-l E 5 



R = 



Aj y |T,(k)| 2 

8 V 



c k 



(29) 



The phase stiffness J is related to the superfluid density n s by 
J = n s /{4M). The superfluid density n s can also be obtained 
from its standard definition ll23ll . When the superfluid moves 
with a uniform velocity v s , the pair field transforms as O — > 
Oe 2lMVl ' T . The superfluid density n s is defined as the response 
of the thermodynamic potential D. to an infinitesimal velocity 
v s ; that is, 0(u s ) = Q(0) + \n s v 2 + 0(v 4 s ). For the present 
system, the superfluid density equals the total fermion density 
n at zero temperature, guaranteed by the Galilean invariance . 

The dispersion of the gapless Anderson-Bogoliubov mode 
is given by 



w(q) = c s |q|, 
where the sound velocity c s reads 



(30) 



n 
Ma' 



(33) 



We note that this relation is quite generic. It also applies to 
s-wave fermionic superfluids. Since n goes smoothly with /j 
we conclude that the sound velocity behaves nonanalytically 
across the quantum phase transition point fi — 0. In Fig. [2] 
we show the behavior of the sound velocity c s around the 
BCS-BEC quantum phase transition, using the same poten- 
tial model as employed in Fig. Q] It shows obviously that the 
evolution of the Anderson-Bogoliubov mode is not smooth, 
corresponding to the BCS-BEC quantum phase transition. 



0.23 




0.22 



0.21 



0.19 



0.18 



0.17 



FIG. 2: The velocity of the Anderson-Bogoliubov mode as a function 
of the density n for the isotropic p x + ip y pairing (upper panel) and the 
anisotropic p x pairing (lower panel). In the calculation we employ 
the same potential model as used in Fig. [JJ All quantities are shown 
in proper units. 



./ 



R + B 2 /A 



(31) 



By comparing the expansion parameters A, B, and R with the 
expressions in Eq. (21), we find that 

1 dn(M, A ) 

K — , 

4 3ju 



B 



1 dn(/j., Ao) 

2 <9A ' 
d 2 V(A ,fi) 

dti ' 



(32) 



III. rf-WAVE PAIRING IN SPIN-i FERMI GASES 

The many -body Hamiltonian of 2D spin-i fermions can be 
written as H — Hq + H mt , where the single-particle part reads 

#0 = Ek&C^k^k + ^[k^ik) an ^ tne d-wave interaction 
part is 

Hint = J] y ^ k ' k 'K,k +q /2 ( /'I.-k +q /2 i /'|.-k' +q /2 i AT,k' +q /2-( 34 ) 
k,k', q 
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The generic infrared behavior of the rf-wave interaction po- 
tential is Vd(k, k') ~ k 2 (k') 2 for k, k! — > 0. Here we also con- 
sider a separable potential for d-wave interaction, V</(k, k') = 
-/l(/r</(k)FJ(k'), where A d is the of-wave coupling constant and 
the function F ( /(k) characterizes the of- wave pairing symmetry. 
The infrared behavior of the function F</(k) is 



r d (k) ~ k 2 , k ^ 0. 



(35) 



As a result, the infrared divergence at // = is more pro- 
nounced for c/-wave pairing. 

The partition function of the system is given by X, — 
j DipDif/'' exp {-J>[^, i/f ' ]}, where the action reads 



l/T 



dT 



k,o- 



■ (36) 



Again, we introduce the auxiliary complex pairing field O q (r), 
which couples to ^^t, and apply the Hubbard-Stratonovich 
transformation. With the help of the Nambu-Gor'kov basis de- 
fined as *F k = (if/ k , tf/'l k ) T we obtain Z = j exp { - 
<S e ff[<E>,<l>*]}, where 



>eff 



\ %(t)\ 2 



^(ft^-TrlnG^,) 



k.k' 



(37) 

Here the inverse single-particle Green's function G k 1 k , takes a 
similar form: 



kk ' 1 ^-k-(T)r*(^) -(d T -&)6 kX 



(38) 



Following the same method used in the previous section, we 
decompose the pairing field as <E> q (r) = Ao^ q ,o + 4>q(r), where 
4> is the fluctuation around the mean field. The effective action 
^ e ff [O, O*] can be expanded in powers of the fluctuation <p, 



<S eff [0, O*] = <S®(Ao)+S^[0, f] + ■ ■ ■ , where S^(A Q ) is the 
saddle point or mean-field effective action with Aq determined 
by the saddle point condition dSf^/dAo = 0. 

A. Quantum phase transition 

The effective potential for rf-wave pairing at zero tempera- 
ture is given by 



J°)/ 



T(Ao,ai)= -r + Y (fk-£k). 



(39) 



Note that the factor 1 /2 in Eq. (10) is replaced by 1 here due to 
the appearance of spin degree of freedom. The single-particle 
excitation spectrum E k is also given by 



£ k = yfil + |A(k)| 2 , 



(40) 



where A(k) = AoF,;(k) is the gap function for rf-wave pair- 
ing. The order parameter Aq is determined by the li-wave gap 



equation 



|F^(k)| 2 
2E k 



(41) 



To study the nature of the quantum phase transition, we also 
focus on the thermodynamic potential Cl(fi) = *y(Ao(/i),//). 
The first derivative of the thermodynamic potential with re- 
spect to n gives the number equation, 



n(p., Aq) = - 



d£l(M) 



Then we consider the second derivative 

d 2 Q.(/u) 



a(p) = -- 



dn 2 



Again, this quantity can be expressed as 

3«Ou,A ) (dn(M,Ao)\ 2 (d 1< V(A ,ti)) 



a(p) 



djj 



dAo 



dA 2 



(42) 



(43) 



(44) 



We note that the second term was missing in the previous 
study lfl6ll . However, the divergent behavior of a is governed 
by the first term. Therefore, the results in Ref. Ill 611 are still re- 
liable. The derivatives in Eq. (44) can be explicitly evaluated 



dnip, A ) 
dfi 

dn(n, A ) 
<3Aq 



ir^(k)| 2 

E 3 ' 
c k 

^k|F rf (k)| 2 



E 3 
c k 



d 2f V(A ,fi) 



M 2 



ir^k)! 4 

E 3 

c k 



(45) 



By using a NSR-type potential, Botelho and Sa de Melo 
showed that a is divergent at = lfl6ll . Here we show 
that this divergence generally appears due to an infrared diver- 
gence at yu = 0. For k — » 0, we have T f /(k) ~ k 2 and therefore 
E k ~ k 2 at fi = 0. As a result, the momentum integral 



IF^k)! 2 

E 3 



(46) 



is infrared divergent at fi = 0. Actually, at = 0, the infrared 
behavior of the above integral is 



J 2 (0) 



dk 

T' 



(47) 



Therefore, the quantity a(/i) itself is divergent at /j — 0, which 
indicates a quantum phase transition. This is not surprising, 
because the infrared divergence should be more pronounced 
for c/-wave pairing. 

The above discussion shows that a quantum phase transi- 
tion occurs at fi = for cZ-wave pairing, indicated by the di- 
vergence of the quantity a(p). The divergence is also caused 
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FIG. 3: The second derivative of the thermodynamic potential Sl(fj) 
with respect to ft as a function of the density n for the isotropic 
d x i_^ + 2id xy pairing (upper panel) and the anisotropic rf v 2_ v 2 pair- 
ing (lower panel). The inserts show the chemical potential fi as a 
function of n. In this calculation, we employ the NSR-type potential 
for d-wave interaction, T(k) = (k/ki) 2 g(<p)/(l + k/k ) V2 [Ql- Here 
g(tp) = cos2i/> for rf v 2_ v : pairing and g(<p) = e 2 '^ for rf r 2_ T 2 + lid^ 
pairing. The density n has been normalized such that n = 1 at fi = 0. 
All quantities are shown in proper units. 



B. Collective modes 

To study the behavior of the collective modes across the 
BCS-BEC quantum phase transition for csf-wave pairing, we 
also consider the Gaussian-fluctuation part S ef[ [<p, (f>*], which 
can also be written as 

5 S = \ X At (e)M(G)A(0- (48) 



The matrix elements of M(0 are also constructed using the 
mean-field fermion propagator @{K) given by 



r~ l (K\ - I ioj " ~ ^ k A ° r '^ k ) 

y (A> -\ A r*(k) /«,„ + & f 



(49) 



We have 

M„(0 = -^- + ^ |rrf(k)|2Tr[ ^ ll(/:+) ^ 22(ii: - )] ' 

d K 

m 12 (0 = ^|r,(k)| 2 Tr[^ I2 (/i: + )^ 12 (/:_)]. (50) 



The relations M 22 (0 = M u (-fl) and M 21 (0 = M 12 (0 
also hold here. By completing the summation over the boson 
Matsubara frequency, we obtain at zero temperature 



M u (<2) 



Ad V \iv n -E + - 



2 2 
V.V_ 



M 12 (0 



E+ - E- iv n 



Xiir^k)! 2 



2 I U+V+U-V- 



iv„ — E+ — E- iv„ +£ , + +£ , _ J 



(51) 



solely by the infrared behavior of the c/-wave pairing potential 
and is independent of the details of the interaction as well as 
the order parameter symmetry (i.e., the angle dependence of 
the gap function). 

In Fig. [3j we show the numerical results of a{ji) for two typ- 
ical cZ-wave pairings: the complex and isotropic d x 2_ x i + 2id xv 
pairing and the anisotropic d x 2_ y i pairing. For d x 2-f + 2id xv 
pairing, the superfluid state is fully gapped for both fi > 
and fi < 0, and it is only gapless at fi = 0. For d x i_ y 2 pair- 
ing, the superfluid state is gapped for fi < and gapless for 
jt/ > 0. In the calculation we have also adopted a NSR-type 
potential to regularize the ultraviolet divergence in the gap 
equation. However, as we mentioned above, the divergence 
at fi = does not depend on this specific choice. The results 
in Fig. [3] show that the divergence appears for both isotropic 
and anisotropic pairings, consistent with our conclusion that 
the divergence is solely due to the infrared behavior of the 
rf-wave pairing interaction. 



Here the notations are the same as in the last section. 

The dispersions cj(q) of the collective modes are also de- 
termined by the equation detM[w(q), q] = 0. Again, we de- 
compose Mn(w, q) as Mn(w, q) = M|j(w, q) + M n (w, q), 
where M|j (a>, q) and M n (a>, q) are even and odd functions 
of a), respectively. Then we separate the complex field (p(Q) 
into its amplitude part A(Q) and phase part 9(Q), 4>{Q) = 
M.Q) + /An 0(0. The effective action «S2! becomes 

^4^(^(0 ^0)N(0(*g), (52) 

where the matrix elements of N(0 read Nn(0 = 2(M[j + 
M 12 ), N 22 (0 = 2A 2 (M[ 1 - M 12 ), N 12 (0 = 2/AqM" , and 
N 2 i(0 = -2/AnM n . At the saddle point we also have pre- 
cisely M|j (0, 0) = Mi 2 (0, 0) for rf-wave pairing. Therefore, 
the phase mode at q = is gapless, corresponding to the 
Anderson-Bogoliubov mode for c/-wave pairing. 

The small q and a> expansion of N(0 also takes the form 
Nn - A + Cq 2 - Do? + ■ ■ • , N 22 = 7q 2 - Rar + ■ ■ • , and 
Ni 2 = -N 2 i = -iBoj+- ■ ■ . The explicit forms of the expansion 
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parameters A, B, D, and R here are given by 

ir d (k)| 4 



E 3 



B = 



D = 



Ap y &|r d (k)p 

k k 



IV 



^ir rf (k)p 



7? 



4 V 

4 V 



ir rf (k)p 

c k 



(53) 



The phase stiffness 7 is also related to the superfiuid density 
n s by J — n s /(4M) with n s = n at zero temperature. 

The dispersion of the gapless Anderson-Bogoliubov mode 
is given by w(q) = c s \q\, where the sound velocity c s reads 



./ 



R + B 2 /A ' 



(54) 



Comparing the expansion parameters A, B, and R with the ex- 
pressions in Eq. (45), we also obtain the following relations 
for of- wave pairing: 



R = - 



1 dn(ji, A ) 
4 



B - — 



1 3n(/i, A ) 
2" 



A = 



dA 



Therefore, for d-wave pairing we also have 



n 
Ma' 



(55) 



(56) 



Since the quantity a is divergent at the quantum phase transi- 
tion point n — 0, the sound velocity c s goes nonanalytically 
across the quantum phase transition and vanishes at ft — 0. In 
Fig. H] we show the behavior of the sound velocity c s around 
the BCS-BEC quantum phase transition, using the same po- 
tential model for of-wave pairing as employed in Fig. [3] It 
shows that the evolution of the Anderson-Bogoliubov mode 
for d wave is also not smooth. The vanishing of the sound ve- 
locity c s may bring interesting thermodynamic consequences. 
For instance, the low-temperature specific heat caused by the 
Goldstone mode can be given by 



2n 2 
15c 3 



(57) 



Therefore, the low-temperature specific heat of the fermionic 
superfluids near ft = should be very large for c/-wave pairing. 

On the other hand, the expansion parameter D, which 
is related to the massive amplitude mode or the so-called 
Anderson-Higgs mode, is divergent at the quantum phase 
transition point fi = 0. A rough estimation of the mass gap of 
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FIG. 4: The velocity of the Anderson-Bogoliubov mode as a function 
of the density n for the isotropic rf,.2_ v 2 + 2id„, pairing (upper panel) 
and the anisotropic d x i_ y i pairing (lower panel). In the calculation 
we employ the same potential model as used in Fig. [3] All quantities 
are shown in proper units. 



the amplitude mode can be given by Mah = + B 2 /R)/D, 
which vanishes at fi = 0. This indicates that the ampli- 
tude mode gets softened around the quantum phase transition 
point, which could be interesting for future identification of 
the massive Anderson-Higgs mode in fermionic superfluids. 



IV. SUMMARY 

While it is generally accepted that the BCS-BEC evolu- 
tion in fermionic systems with s-wave pairing is a smooth 
crossover, for nonzero orbital-angular-momentum pairing 
such as p- or of-wave pairing, the BCS-BEC evolution is asso- 
ciated with a quantum phase transition at vanishing chemical 
potential yu. In this paper, we have studied some generic fea- 
tures of the BCS-BEC quantum phase transition and the col- 
lective excitations in two-dimensional fermionic systems with 
p- and cZ-wave pairings. Our generic conclusions can be sum- 
marized as follows: 

(1) The quantum phase transition in two-dimensional 
fermionic systems is essentially related to the infrared behav- 
ior of the pairing interaction, which causes the infrared be- 
havior of the fermionic excitation at fi = 0: ~ k 1 , where 
I = 1 for p-wave pairing and I = 2 for rf-wave pairing. The 
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nonanalyticities of the thermodynamic quantities are due to 
the infrared divergence caused by the fermionic excitation at 
H = 0. 

(2) The evolution of the Anderson-Bogoliubov mode is not 
smooth: Its velocity is nonanalytical across the quantum 
phase transition, due to the infrared divergence caused by the 
fermionic excitation at /i = 0. 

(3) The BCS-BEC quantum phase transition and nonsmooth 
evolution of the collective modes in 2D systems with nonzero 
orbital-angular-momentum pairing are solely related to the in- 
frared behavior of the pairing interaction and are independent 
of the details of the interaction potential as well as the pairing 



symmetry. 

Finally, we point out that while we have studied the 2D con- 
tinuum models, the generic features of the BCS-BEC quantum 
phase transition summarized above also apply to fermions in 
a 2D square lattice 11121 12211 since the infrared behavior of the 
fermionic quasiparticles remains. Therefore, it is interesting 
to extend our studies to lattice systems. 
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